Quantum phase transition in a three-level atom-molecule system 
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We adopt a three-level bosonic model to investigate the quantum phase transition in an ultra- 
cold atom-molecule conversion system which includes one atomic mode and two molecular modes. 
Through thoroughly exploring the properties of energy level structure, fidelity, and adiabatical 
geometric phase, we confirm that the system exists a second-order phase transition from an atom- 
molecule mixture phase to a pure molecule phase. We give the explicit expression of the critical 
point and obtain two scaling laws to characterize this transition. In particular we find that both the 
critical exponents and the behaviors of ground-state geometric phase change obviously in contrast to 
a similar two-level model. Our analytical calculations show that the ground-state geometric phase 
jumps from zero to 7r/3 at the critical point. This discontinuous behavior has been checked by 
numerical simulations and it can be used to identify the phase transition in the system. 

PACS numbers: 05.30.Rt, 67.85. Hj, 03.65.Vf 



I. INTRODUCTION 

Quantum phase transition (QPT) is one of the most 
important concepts in the many-body quantum theory. 
As a central fundamental transition phenomenon at the 
temperature of absolute zero, it describes an abrupt 
change in the ground state of a many-body system due 
to its quantum fluctuations P, . The experimental ob- 
servation of a QPT from a superfluid (SF) to a Mott 
insulator (MI) in a gas of ultracold atoms [sl] inspired 
great interest in investigating clean, highly controllable, 
and strong correlated bosonic systems [1]. Actually, the 
ultracold atomic gases P] have become an ideal platform 
to study many-body physics because of their enormous 
applications and the advanced experimental techniques 
available in the fields of atomic and optical physics [6|. 

In recent years, a remarkable development in the afore- 
mentioned filed is to convert ultracold atoms to molecules 
via Feshbach resonance 0-[ll or photoassociation (lol - 
[T^ techniques. Compared with the fermionic model in 
this kind of systems, the bosonic model is of interest for 
theoretically exploring the QPTs. When both ultracold 
atoms and molecules are bosons, the systems possess a 
few degrees of freedom and a large particle number which 
can greatly simplify the calculations. On the other hand, 
the atom-molecule conversion systems can be well de- 
scribed by a mean-field theory when the particle number 
is large enough and this treatment will lead to nonlin- 
earity. These features have stimulated much efforts to 
study the adiabatic evolutio n [isl-fisl . geometric phase 
[l^ll^l, and phase transition 21, 22l| of the systems. In- 
stead of the traditional approaches for describing a QTP 
(i.e., using the concepts of order parameter and symme- 
try breaking), very recently Santos el al. adopted the 
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concepts of entanglement and fidelity to investigate the 
QPT in a two-level bosonic atom- molecule system j2l| . 
Motivated by this work we discussed same problems from 
the perspectives of scaling laws and Berry curvature . 
However, the connection between the mean- field Berry 
phase and the phase transition in this type of bosonic 
model, and the properties of the QPT in a three-level 
atom-molecule system are still unresolved, which call for 
further theoretical considerations. 



As a continuous work, in this paper we investigate 
the quantum phase transition in an ultracold atom- 
molecule conversion system by adopting a A-type three- 
level bosonic model. Based on this model, we first discuss 
the structure of quantum energy levels and analyze the 
properties of ground state. In order to compare the phase 
transition properties with a similar two- level model [2]| , 
we study the energy gap, the ground-state fidelity, and 
mean-field geometric phase. We illustrate that, when 
the ratio of the coupling strength between two molecu- 
lar modes to that between atomic mode and the upper 
molecular mode exceeds a critical value, a similar QPT 
from a mixed atom-molecule phase to a pure molecular 
phase is also observed in our system. To characterize this 
transition we obtain the analytical expression of the crit- 
ical point by using the mean-field approach and derive 
two critical exponents via numerically studying scaling 
laws. In particular we calculate the ground-state geo- 
metric phase and find its discontinuous behavior at the 
phase transition point. 



Our paper is organized as follow: In Sec. II we give the 
second-quantized model and its mean-field description. 
In Sec. Ill we explore the properties of energy levels and 
ground states. In Sec. IV we choose the characteristic 
scaling law, fidelity, and adiabatic geometric phase to 
describe the QPT. Section V presents our conclusion. 
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II. THREE-LEVEL MODEL AND MEAN-FIELD 
DESCRIPTION 

The system we consider here is illustrated schemati- 
cally in Fig. [T] It describes the process of creating ultra- 
cold diatomic molecules from bosonic condensed atoms, 
which constitutes a A-type three-level model. In the 
three- mode description, each mode \a) {a — a,g, and e 
respectively represent the atomic mode, the ground-state 
molecular mode, and the excited-state molecular mode) 
is associated with an annihilation operator /? (/? = a, 6g, 
and be) due to the basic assumption that the spatial 
wavefunctions for these modes are fixed. By setting the 
energy of atomic mode as zero, the Hamiltonian of the 
system takes the following second-quantized form with 
S= 1 I20I: 



-blaa + R.c, (1) 



N 

where the abbreviation H.c. denotes the operation of Her- 
mitian conjugate. Vd and Vp are the frequencies of two 
laser pulses rid and Qp, respectively. The frequencies 
ujg and We measure the molecular ground state and ex- 
cited state energies, respectively. The total atom number 
N ^ Na + 2{Ng + Ne) with Na = flta, Ng = blbg, and 
Nf, = bl,be, commutes with the Hamiltonian ([1} and is 
therefore conserved. Notice that the laser pulse param- 
eter rip can be complex. To achieve this one can split 
the laser pulse into two beams and then recombine and 
focus them on the system. As a result, we can express 
the complex parameter flp as rip = + ^2^~^^ with 1^1 
and ^2 being real numbers, where the phase factor ip is 
determined by the difference of optical paths between the 
two laser beams [20| . 

For convenience, we rewrite the above Schrodinger pic- 
ture Hamiltonian as, Hs = Hq + Hi, where 



Ho =Vpb\be + {vp - Vd)'blbg, 

Hi ={uje ~ Vp)blbe + {ojg -iyp + yd)blbg+ 



ride^^^'blL + 



N 



-blaa + H.c. 



(3) 



then we choose uje 



Vd and apply the interaction 



picture, i.e.. Hi — e'^^°*Hie *^ot^ i\yQ Hamiltonian fi- 
nally becomes 




FIG. 1: Schematic diagram of a three-level atom-molecule 
conversion model coupled by two laser pulses Qd and Sip. 
The stable atomic state, the molecular ground state, and the 
molecular excited state are denoted by |a), |g), and |e), re- 
spectively. A measures the detuning of the pump laser fip 
with respect to the transition from |a) to |e). 



proven to be a powerful tool for studying ultracold atoms 
and Bose-Einstein condensates (BECs). In the semiclas- 
sical limit N ^ 00, the quantum model becomes clas- 
sical and one can replace the operator 13 with a corre- 
sponding complex number /3 (/3 = a,bg,be), i.e., "H = 
A{\be\^ + \bg\^) + z{b:bg + b;be) + p[e'''^b:a^ + e'He{a*)^]. 
By using the equations idp/dt = dH/dp*, we can ob- 
tain the following Schrodinger equations together with 
the normalization condition |ap + 2(|6gp + l^eP) = 1 to 
govern the dynamical behaviors of the system. 



where 



H 



mf 









(5) 



(6) 



(■2) and 1-0) = {a,bg,be)"^ ■ It is worth emphasizing that. 



although the collisions between ultracold particles have 
been neglected in our model, the nonlinearity also arises 
from the mean-field treatment for the fact of two atoms 
to form one molecule. Mathematically, the mean-field 
Hamiltonian Hmf is a function of the instantaneous wave 
function as well as its conjugate. It is a nonhermitian ma- 
trix and it is invariant under the following transformation 



Hi = A{b% + blbg) + zblbg + ^^blaa + H.c, (4) 



-icf) ^ 



N 



e 
26 
29 



IV'). (7) 



where the new parameters A 



UJe 



and = aTg{Qp) have been introduced. 

To complement the quantum description and gain in- 
sights into the existence of a QPT in our model, we 
adopt a semiclassical description of the system by fol- 
lowing the usual mean-field approach, which has been 



The lack of U{1) gauge transformation is a particular 
interesting point of the above mean-field model, which 
may lead to some new properties of the system. In the 
subsequent sections, based on models (H)) and we 
will discuss the QPT in the system both from the fully 
quantum perspective and the mean-field perspective. 
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z z versus the parameter z with A = and A = 0.5. Both z and 

A are rescaled by p. 

FIG. 2: (Color online) Quantum energy levels for different 
particle numbers: (a) = 2, (b) iV = 4, (c) = 6, and (d) 
N = S. The thick and thin lines denote the cases A = and 
A = 0.2, respectively. The energies are shown in units of Np 
while the parameters z and A are rescaled by p. 



III. ENERGY LEVELS AND GROUND STATES 

Taking advantage of the fact of N is conserved, one can 
diagonalize the quantum Hamiltonian Hj. For simplicity, 
hereafter we assume that N takes an even constant value, 
then the Hilbert space of the A^-particle system can re- 
duce to + + 2) dimension in the Fock basis, 

i.e., \na)\ng)\ne) = (na!n5!ne!)-i/2(at)"<. (Sj |0) 
with |0) being the vacuum state, where Ug — 0,1, ■■ ■ , 
Tie = 0, 1, • • • , Y ~ fig, and Ua ^ N — 2{ng + n^) repre- 
sent the populations of particle in states \bg), \be), and 
I a), respectively. By directly diagonalizing the Hamilto- 
nian matrix with a fixed N, we obtain the eigen-energy 
levels and the ground states of the system as shown in 
Figs. [Hand 131 respectively. 

It is well known that a typical A-type three-level sys- 
tem supports dark-state solutions with zero eigenvalue 
[23l - [25j . This type of state can result in a phenomenon 
known as coherent population trapping (CPT). For our 
system, when A = 0, from Figs. [Hb)-md) we see that the 
energy levels with energy value being zero are degenerate 
while other nonzero-energy levels are nondegenerate and 
are symmetrically distributed in both sides of the center 
level. This symmetrical energy structure is determined 
by the symmetry of the Hamiltonian Hj with (f> = 0, i.e., 
the change of variables {z,p) — >■ —{z,p) is equivalent to 
the unitary transformation be — >■ —be- In this case, the 
degeneracy of the zero-energy level (i.e., d) is given by 

d= r(^ + l)/2l =i(7V-Mod[7V,4]) + l, (8) 



the symbol [ ] stands for the ceiling function which maps 
a real number to the smallest following integer. 



Notice that, if the parameter A has a perturbation, the 
above mentioned symmetry of the system with A = will 
be broken. This leads to the energy levels shift and the 
zero-energy level splitting. For example, when A = 0.2 
and N — 8 [see Fig. [TJd)], all energy levels have been 
pushed up and the zero-energy level has split into three 
nonzero-energy levels. For different N, the maximum 
number of energy levels should be (y + 1)(y + 2)/2- 



Now we discuss the ground-state properties which are 
closely associated with the QPT in the system. On one 
hand, by diagonalizing the Hamiltonian Hi numerically, 
both for A = and for A 7^ 0, we have calculated the 
ground states with different total atom numbers. The 
results for atomic population fraction (i.e., Na/N) in 
the ground state are demonstrated in Fig. |31 We find 
that the atomic fraction in the ground state decreases 
and gradually approaches to zero as the ratio of the 
coupling strength between two molecular modes to that 
between atom mode and the molecular mode increases. 
On the other hand, we study the ground state from 
the mean-field perspective. Based on the model (O, we 
solve the mean-field ground state from the eigen-equation 
Hmf{a*,a)\'ip) = 9(/x)|'(/') with /i being the chemical po- 
tential for atoms and l-tp) being the eigenstate. For A = 0, 
z > 0, and p > 0, we obtain the eigenvalue and the cor- 
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state (a) and its first (b) and second (c) derivatives with re- 
spect to coupling z. The parameters z and A are rescaled by 

p FIG. 5: (Color online) (a) The total atomic number A'^ versus 

the offset between zn and Zc- (b) The minimum value of the 
energy gap (i.e., AiSmin) at the pseudo-critical point zn versus 
A^. Solid lines are plotted for guiding eyes, respectively, (from 

responding eigenfunction for the ground state as follows: left to right) with the slopes being -1.54764 and -0.32912. 




, Z < 2/9. 



(9) IV. QUANTUM PHASE TRANSITION 

The previous calculations and analysis have demon- 
strated that the process of converting ultracold atoms 
to homonuclear diatomic molecules in bosonic system 
(]^Q) is a QPT which differs from the well-known BCS-BEC 
crossover phenomena in fermionic systems (26l . [13] . Sub- 
sequently, we will describe and characterize this phase 
transition from different perspectives. 



For A 7^ 0, although the solutions to ground state can 
also be obtained analytically, the expressions are gener- 
ally too messy to be instructive. We therefore simply dis- 
play the results in Fig. [31 From Fig. [3]we find that, both 
for A = and for A ^ 0, the results for the quantum 
model (|H) will tend to the analytical mean-field results 
with increasing the total particle number N . 

It must be mentioned that, for our nonlinear system 
([5]) , the classical energy £ does not equal the chemical po- 
tential, and the relation between them is £" = /x±p|6e| |ap. 
We have calculated the ground-state energy £o analyti- 
cally and find its second derivative possesses a disconti- 
nuity at a critical point Zc = 2p -f A as demonstrated in 
Fig. m This divergence behavior implies that the system 
exists a second-order phase transition in the thermody- 
namic limit. When z < Zc the system is in an atom- 
molecule mixture phase (i.e., |ap > 0) and when z > z^ 
the system is in a pure molecule phase where |ap = 0. In 
the mixture phase, the asymptotic behavior of the ground 
state in the vicinity of the critical point with A = is 
given by the variation of the parameter sq = |ap, i.e., 

soU^Ze = ^[4-Zc(2z-Zc)]. (11) 



A. Scaling laws 

In order to understand the QPT in the system, we be- 
gin our discussion by analyzing the dimensionless energy 
gap between the first excited state and the ground state, 
namely, /S.E = {Ei — Eq)/ p- With the help of diagonaliz- 
ing the Hamiltonian ^ numerically, we have calculated 
the energy levels with different particle numbers. For a 
fixed N, the energy gap takes a minimum value at a point 
zat (it can be viewed as a pseudo-critical point of the N- 
particle system) and this point just corresponds to the 
position of avoided level crossing (see Fig. [5]). Generally, 
the QPTs often occur at the positions of level crossings 
or avoided level crossings. For our avoided level crossings 
system, the existence of the minimum of energy gap indi- 
cates a basic signature of the phase transition. Similar to 
the phenomena studied in a two-level atom-molecule sys- 
tem [2l| and in other systems, we find that, the gap AE 
in our system also tends to zero at a single point rather 
than over an interval of the dimensionless parameter z/p 
with the particle number N ^ oo. This specific phe- 
nomenon implies that, when N ^ oo, the ground state 
is degenerate at the point Zc where is no phase, which is 
a requirement for the occurrence of a broken-symmetry 
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phase (2l|. 

To capture more features of the QPT in the system, 
we study the scaHng behavior of the energy gap near the 
critical point. To this end, we first calculate the energy 
gap for different parameter A and the results have been 
plotted in Fig. [S] Either for the variation of the to- 
tal atomic number N versus Iz^r — [see Fig. [D^a)] or 
for the change of the minimum of the energy gap (i.e., 
Ai^min) with respect to [see Fig. [Sfb)], the same char- 
acteristic scaling laws have been observed for different A, 
and different lines in each figures with a same slope gives 
the evidence. 

In the quantum model (j4]), the total atom number N 
can be regard as a correlation length scale of the system, 
and then one can connects this length scale to the offset 
between the pseudo-critical point and the critical point. 
Quantificationally, we have 



zn\ 



N- 



(12) 



where 1/ ~ 1.54764 is a critical exponent and k ~ 0.18273 
is a inessential constant. This scaling law shows that the 
pseudo-critical point changes and tends as N~^^'^ toward 
the critical point and clearly approaches N 00 

[see Fig. [SJ^a)]. From Fig. ^h), we find another scaling 
law, that is 



(13) 



where F ~ 1.67506 is a constant. C ~ 1.32912 gives an- 
other important exponent, namely, the dynamic critical 
exponent. It must be mentioned that all constants and 
exponents given in the above two formulas are obtained 
in the case of A = 0, for other cases their values may 
have a slightly change. Comparing the product of two 
exponents in our swtem with that in a two-level atom- 
molecule system [23 we find that the values are obvi- 
ously different. This difference indicates that the above 
two models belong to different universality classes. 



B. Fidelity 

Similar to other concepts, the behavior of the fidelity 
can also be employed to identify the phase transition 
[H, [2^. The fidelity is a measure of the distance be- 
tween two states and this concept has been widely used 
in the field of quantum information [30| . One can define 
the fidelity through the modulus of the wave functions 
overlap between two states, i.e.. 



F(,Al,V'2) = |(V^l|^2)|. 



(14) 



Here we only focus on the behaviors of the fidelity be- 
tween two ground states. One ground state is obtained 
when A = 0, namely, |A = 0), the other ground state 
is calculated by treating A as a perturbation parameter, 
denoted as |A ^ 0). we have estimated the wave-function 
overlap between two ground states with different parti- 
cle number N and varying parameter A. Figure [5] shows 




FIG. 6: (Color online) Ground-state fidelity versus the pa- 
rameter z with A*' = 50 (a) and with A'^ = 100 (b) for diflferent 
A. The parameters A and p are rescaled by p. 



the fidelities between the ground state |A = 0) and the 
ground state |A — a) with a — 0.1,0.2,0.3, and 0.4. 
Both for iV = 50 and for N = 100 we see that the fidelity 
|(A — 0|A 7^ 0)1 shows dip behavior at the point corre- 
sponding to the pseudo-critical point. The results imply 
that the two ground states are distinguishable and there 
is an obvious signal for the QPT as long as a ^ 0. More- 
over, we observe that the dip of the ground-state fidelity 
becomes deeper and the point where the fidelity has a 
minimum varies with increasing the value of a, this phe- 
nomenon is very different from that in a two-level atom- 
molecule model where the fidelity has a minimum at the 
same point [2l| . The reason is that, in our model the 
phase transition point Zc is a function of the parameters 
p and A. For a larger value of a, we find that the distin- 
guishability of the two states increases and the minimum 
of the fidelity changes evidently. 



Now we compare the results obtained in the case of 
TV = 50 with the results for N — 100. For a same A, 
it is seen that the position of the minimum fidelity for 
N ~ 100 is closer to the critical point Zc than that for 
= 50. In fact we have calculated the ground-state fi- 
delity with varying the particle number N and find that, 
for a fixed a, with increasing A^, the wave functions over- 
lap between the two states become smaller (i.e., the two 
states are more distinguishable) and the position where 
the minimum fidelity occurs moves toward 2p + A. Thus 
in the finite particle number case the occurrence of the 
minima of the fidelity gives a information about the phase 
transition in the system. 
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C. Geometric phase 

In this subsection we convert to investigate the behav- 
ior of the ground-state geometric phase starting from the 
mean- field model ([S]). In the following study we only 
focus our attention on the situation that the detuning is 
absent (i.e., A = 0). Because the analytical results can be 
obtained in this case. To employ the procedures for calcu- 
lating geometric phase in nonlinear systems proposed in 
Refs. [1^ HH, we introduce some new variables through 
setting a = ^1 - 2(pi -f p2)e*^, &g = y^e^^^^+'i', and 
be = -y/pje''-^'*'^'^'', where A = arg(a) denotes the to- 
tal phase, pi = and p2 — |&eP are the popula- 



tion probabilities of the ground state and excited state 
molecules, respectively, qi = arg(6g) — 2arg(a) and q2 = 
arg(6e) — 2arg(a) measure the relative phases. With the 
help of these new variables, the three-level system can be 
cast into a classical Hamiltonian 



n ^2zy/pip2 cos(9i - q2) 

+ 2/oVP2[1 - 2{pi + P2)] cos(g2 



(15) 



The Schrodinger equations ([5]) together with the normal- 
ization condition lead to 



dX 



'2pVP2Cos(g2 + 0), 



and 



dpi _ dV. dqi _ dTi 
dt dqi ' dt dpi ' 



(16) 



(17) 



with i — 1,2. These four equations have estab- 
lished a connection between the projected Hilbert space 
spanned by S(pi, qi) and the parameter space spanned by 
R(2;,p, (/)). 

Now we calculate the geometric phase for the ground 
state of the system. For simplicity, we construct a closed 
loop C in the parameter space by treating z and p as con- 
stants and varying (j) with time from to 27r. The system 
is assumed to evolve adiabatically along the cyclic path 
with a rate e ^ |^| ~ i, where T is the period of the 
cyclic evolution. Initially, we prepare the system in the 
ground state of H„if{(f) = 0), after a cyclic adiabatic evo- 
lution, the state will acquire a geometric phase besides 
the dynamical phase. Here it is noted that the adia- 
batic parameter e 1, and thus it can be regarded as a 
small parameter during the process for determining the 
geometric phase. Following the method in Ref. [3l|, we 
expand the total phase A in a perturbation series in e. 



^ = Ao(60) + Ai(6i) + 0(62), (18) 

to separate the pure geometric part from the total phase. 
The time integrals of the zero-order term and the first- 
order term in Eq. (jl8p will respectively give the dynamic 
phase and the geometric phase in the adiabatic limit e — >■ 
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FIG. 7: (Color online) Comparison between the numerical 
results for the ground-state geometric phase with T — 500 
and the analytical results given by Eq. (|23[l . The dotted hne 
denotes the results obtained from Berry's formula. The inset 
illustrates the numerical results around the critical point with 
different T. The parameter z is rescaled by p. 



or T — 00, and the contribution from the higher-order 
terms will vanish. 

During the adiabatic evolution the system will fiuc- 
tuate around the ground state due to the small but fi- 
nite value of e. This fact allows us to make expansions 
as Pi — Pi(R) + (5pi(R) and qi = gi(R) + SqiCR) with 

1 ~ 1,2, where {pi, qi) stand for the instantaneous ground 
state, Spi and Sqi denote the fiuctuations induced by the 
slow change of the system. Substituting these expressions 
back into Eq. ()16|) . when z < 2p, we have 



Ao 



/^o(R-), 



(19) 



^ -2p2y/p^Z + p[p2{6p2 + 6pi -l)+ dp2] ^^Q) 



V-P2 

where the chemical potential of the ground state can be 
expressed as ^0 = —2z^pip2 — 3p^/p2[l — 2{pi + p2)]- 
Moreover, from Eqs. P7|) . we have 



5p2 



2y/p^{zp2 + zpi - Appl^"^) d(j) 
p[z{l + 6p2 + 6pi) — lQpp\^'^] 



(21) 



To deduce Eqs. (PO]) and (PT|). we have used the fixed- 
point equations f|^|(pi.qi) — and the condition -^Sqi ^ 
0{e^). Combining Eq. ^ with Eq. and then 

using the fixed-point values corresponding to the ground 



state, i.e., (pi, gi;P2, 52) = ( 
we obtain 



16p' 



■,g2±7r: 



±7r- 



ld(j) 
6~dt' 



By integrating Eq. 



Ai - 9^- (22) 
over T with respect to time, we 
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FIG. 8; (Color online) The convergence of population in the 
ground-state molecular mode (a) and the geometric phase (b) 
with respect to the evolution period T iox z/ p = 1. 



have 



- I Xidt - 



TT 

3' 



(23) 



As a comparison we give the resuh that is directly calcu- 
lated from the Berry's formula |32[, i.e., 



(24) 



Now we consider the situation when z > 2p. In this 
case, the ground state, i.e., pi = P2 = 1/4, is independent 
of the parameter R. Simple calculations from Eqs. (fTC)) 
and P7)) lead to 



dX z 

= 2 = 



(25) 



of the ground state jumps from zero to tt/B when the sys- 
tem undergoes the phase transition from a mixture phase 
to a pure molecular phase. These analytical predictions 
have been checked by numerically solving Eq. ([5|) or Eqs. 
p6| and ([HI) as illustrated in Fig. [71 We see that, if the 
evolution period T is large enough, the simulated results 
show a good agreement with the analytical predictions. 
We use Fig. |S] to exhibit the convergence behaviors of 
the ground-state and the geometric phase with increas- 
ing T, and a large convergence rate has been observed. 
It is worth emphasizing that, the different values of the 
geometric phase in different parameter regions can be an 
evident signature of the QPT in the system. 

V. CONCLUSION 

We have investigated the quantum phase transition in 
a three-level atom-molecule conversion system. By us- 
ing different approaches we show that the system exists 
a second-order phase transition similar to a QPT exhib- 
ited in a two-level bosonic model [2l|- Firstly, though 
analyzing the properties of energy gap we derive two 
scaling laws and the corresponding critical exponents. 
It is found that the two-level model and our model be- 
long to different universality classes. Secondly, we dis- 
cuss the ground-state fidelity. A minimum value of the 
fidelity near the critical point has been found. Finally, 
we have calculated the ground-state geometric phase and 
a discontinues behavior at the critical point has been 
observed. This phenomenon is similar to that studied 
in a system of a Bose-Einstein condensate in an opti- 
cal cavity [sl]. It establishes a connection between the 
ground-state geometric phase and the QPT in an inter- 
acting atom-molecule bosonic model following the early 
works in spin-chain systems (s^ . Issj . In summary, we 
have demonstrated novel characteristic scaling laws and 
abrupt change of the ground-state geometric phase in the 
vicinity of the critical point, which give the pronounced 
signals toward the existence of a QPT in the system. 



and then 
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